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Overlap fermions implement exact chiral symmetry on the lattice and are thns an appropriate tool for investi¬ 
gating the chiral and topological structure of the QCD vacuum. We study various chiral and topological aspects 
on Liischer-Weisz-type quenched gauge held conhgurations using overlap fermions as a probe. Particular emphasis 
is placed upon the analysis of the spectral density and the localisation properties of the eigenmodes as well as on 
the local structure of topological charge huctuations. 


1. Introduction 

The understanding of the vacuum structure of 
QCD is an important goal of elementary par¬ 
ticle physics. The distribution of topological 
charge density is believed to describe nonpertur- 
bative phenomena like the large rj' mass, the axial 
U(l) anomaly, the 6 dependence and the sponta¬ 
neous breaking of chiral symmetry. In the tra¬ 
ditional instanton picture the excitations of the 
QCD vacuum are phenomenologically modelled 
as an interacting ensemble of instantons and anti- 
instantons. 

Recent advances in implementing chiral sym¬ 
metry on the lattice have made it possible to di¬ 
rectly clarify from first principles to what extent 
the instanton picture of localised solutions of the 
Euclidean Yang-Mills equations of motion realis¬ 
tically describes the QCD vacuum. 

Overlap fermions m have an exact chiral 
symmetry on the lattice 0 and are the cleanest 
known theoretical description of lattice fermions. 
Their implementation of chiral symmetry and the 
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possibility to exactly define the index theorem on 
the lattice allow us to investigate the relation¬ 
ship of topological properties of gauge fields and 
the dynamics of fermions. A further advantage of 
overlap fermions, in contrast to Wilson fermions, 
is that they are automatically 0{a) improved ^ 
and not plagued by exceptional configurations. 

The overlap operator is defined by 

where we use the Wilson-Dirac operator Dw as 
the kernel of the overlap operator. 

The operator has n_ -I- n_|_ exact zero 
modes, DNipn = 0, n_ (n+) being the num¬ 
ber of modes with negative (positive) chirality, 
Tsf/’n = -V'n (75f/’n = +ipn) ■ The index of Dn 
is given hy Q = n_ — n+. The nonzero modes 
A with Djsfipx = \ip\ come in complex conjugate 
pairs A and A* and satisfy = 0. 

The mass parameter p is chosen to be 1.4, 
which is a compromise between the physical re¬ 
quirement of good locality properties of the over¬ 
lap operator and the performance requirement of 
a small condition number of Dw- 
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To compute the ‘sign function’ 


sgn(X) = = 75 sgn(i7), H = 75 X, (2) 

we use the minmax polynomial approximation . 
We compute the 0(10) lowest eigenvalues of H 
and treat the sign function on the corresponding 
subspace exactly. 

Topological studies using the Wilson gauge 
field action suffer from dislocations and should 
be treated with caution. We therefore use the 
Liischer-Weisz gauge field action |3, which sup¬ 
presses dislocations and greatly reduces the num¬ 
ber of unphysical zero modes. 

The Liischer-Weisz gauge field action is given 
by 
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with coefficients ci, C 2 (cq -I- 8 ci -I- 8 c 2 = 1) taken 
from tadpole improved perturbation theory 

To investigate the volume dependence of our 
data we have simulated at 3 different volumes at 
fixed coupling /3 = 8.45. To explore the a depen¬ 
dence of the results, we also employed a 12^ x 24 
lattice at /3 = 8.10 with approximately the same 
physical volume as the 16^ x 32 lattice. 

In the following table we list the lattices used 
in our calculations, where we take rp =0.5 fm to 
set the scale: 


p 

a [fm] 

V 

confs 

modes 

8.45 

0.095 

12^ X 24 

500 

G(50) 

8.45 

0.095 

16^ X 32 

267 

G(140) 

8.45 

0.095 

243 X 48 

186 

G(160) 

8.10 

0.125 

123 X 24 

254 

0(140) 


2. Density and locality of the eigenmodes 
of the overlap operator 

The spontaneous breaking of chiral symmetry 
by the dynamical creation of a nonvanishing chi¬ 
ral condensate, (^^t), is related to the spectral 
density p{X) of the Dirac operator near zero by 
the Banks-Casher relation (Hj (^^t) = —Tr/Vp{0). 

The eigenvalues of lie on a circle around 
(/o,0) with radius p in the complex plane. In 
our spectral analysis we use the improved oper¬ 
ator = (I - aDN/2p)-^DN rmrm , which 
projects the eigenvalues of stereographically 
onto the imaginary axis. The nonzero eigenvalues 
of come in pairs ±iA, while the zero modes 
are untouched. 

Using the Arnoldi-algorithm, we compute the 
0(50) — 0(160) lowest eigenvalues on every con¬ 
figuration. 

The spectral density of the ‘continuous’ modes 
is formally given by 

= ( 4 ) 

A 

where the sum extends over the nonzero eigenval¬ 
ues ±iA of 

Fig.n shows the spectral density together with 
a fit^ using the prediction of chiral perturbation 
theory. For small eigenvalues one can see a strong 
volume dependence of the spectral density. 

In a finite volume, the spectral density is given 
in quenched chiral perturbation theory imsi as 

p{X) = Eeff(U, A) ^ u;(Q)pQ(USeff(U, A)A), (5) 

Q 

where pq{x) = f (J|2 q| (x) - J|q|+i(x) J|Q|_i(a:)) is 
the microscopic spectral density m in the sector 
of fixed topological charge Q defined in terms of 
Bessel functions Jn{x). 

The effective chiral condensate is given by 
Eeff(U,A) = S(1 - 167r25GV(2EA//2)), where <5 
is the chiral log parameter, / denotes the pion 
decay constant and Gy(M^) = dM^GviM"^) 
is the derivative of the pseudo-Goldstone bo¬ 
son propagator without zero momentum modes 

^To consider the effects of higher orders in chiral pertur¬ 
bation theory we added a term aiA -|- a 2 \^ to the fitted 
formula 0. 
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Gv{M^) = 

Taking the weight w{Q) of the sector of topo¬ 
logical charge Q with w{Q) = 1 from our 
‘measured’ charge distributions, we obtain 5 = 
0.20(3) and E = (212(18) MeV)^ (E^(2GeV) = 
(231(20) Me V)3). 



Figure 1. The spectral density at /3 = 8.45 to¬ 
gether with a fit using the finite volume predic¬ 
tions of quenched chiral perturbation theory. 


A useful measure to quantify the localisation of 
eigenmodes | 15lltilT7| is the inverse participation 
ratio (IPR) 


I = vY^p{x)\ 

X 


( 6 ) 


with the scalar density p(x) = ilj\{x)ipx{x) for 
normalised eigenfunctions J2xPi^) — 1- While 
I = V if the density has support only on one 
lattice point, I decreases as the density becomes 
more delocalised, reaching 7=1 when the density 
is maximally spread over all lattice sites. 

Fig. 12 (a) and (b) show the IPR for all eigen¬ 
modes on all considered configurations as a func¬ 
tion of the imaginary part of the eigenvalue on the 
12^ X 24 and 16^ x 32 lattice at /3 = 8.45. One can 
see that the zero modes and the first nonvanishing 
modes are, especially on the large lattice, highly 
localised, while the bulk of the higher modes is 
nonlocalised with 7 < 2. In Fig. [21(c) we plot the 


IPR averaged over bins with bin size AlmA = 50 
MeV, except for the zero modes which are con¬ 
sidered separately. 

From the volume dependence of the IPR one 
can derive the effective dimension d of the 
eigenmodes based on the relation 7 oc 
Fig. m (d) shows the volume dependence of the 
averaged IPR at fixed lattice spacing. We fit 
the zero modes and lowest nonzero modes with 
0 < |Im A| < 50 MeV separately and find d= 
2.3(1) and d=3.0(2), respectively. 

While the low modes strongly depend on V and 
a and are thus localised with an effective dimen¬ 
sion d = 2... 3, the bulk of the higher modes is 
independent on V and a, and extends in d = 4 
dimensions. 

3. The local structure of topological charge 
fluctuations 

The topological charge density for any 
75 -Hermitean Dirac operator satisfying the 
Ginsparg-Wilson relation is defined as m 

dix) = ^TT:"f5D{x,x), Q = ^q{x). (7) 

X 

To compute the topological charge density, we 
use two different approaches nwM . In the first 
approach, we directly calculate the trace of the 
overlap operator according to Eq. Q. This is 
a computationally very demanding task and is 
therefore performed on only 53 (5) configurations 
on the 12^ x 24 (16^ x 32) lattice so far. 

The full density computed in this way includes 
charge fluctuations at all scales. 

The second technique involves the computation 
of the topological charge density based on the low 
lying modes of the overlap Dirac operator. This 
approach is gauge invariant and leaves the lat¬ 
tice scale unchanged, in contrast to the cooling 
method. Using the spectral representation of the 
Dirac operator, the truncated eigenmode expan¬ 
sion of the topological charge density reads 

dXcntix) = - ( 8 ) 

I A| Acut 

C^(x) = lp{{x)^5lPx{x) . 
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(a) 




(b) 



(c) 


V |£m*| 
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Figure 2. The IPR of overlap eigenmodes (for all computed eigenmodes on all considered configurations) 
on the 12^ x 24 lattice (a) and the 24^ x 48 lattice (b) at /3 = 8.45. The averaged IPR is shown in (c), 
while the volume dependence with fixed (3 = 8.45 for the zero and the low lying modes is shown in (d). 


Truncating the expansion acts as an ultravio¬ 
let filter by removing the short-distance fluctua¬ 
tions from q{x). Note that the topological charge 
Q = 9Acut (^) is not affected by the level of 
truncation. 

Let us first consider the local structure of the 
topological charge fluctuations, by studying sign 
coherent clusters formed by connected neighbours 
Xi with \q{x^)\ > qcut = (0.1... 0.4) gmax, where 
gcut is a cutoff value to be varied at will. To 
give a picture of the spatial distribution of the 
topological charge density, isosurface plots with 
q{x) = ±gcut are shown in Fig. O It is obvious 
that the UV cutoff in the topological density re¬ 
sults in a different topological structure. This is 
the structure actually seen by the lowest modes. 

Reflection positivity, i.e. positivity of the met¬ 
rics in Hilbert space, demands that the topologi¬ 
cal charge correlator is negative EH 


Cq{r) = ^'^{q{x)q{y)) <0, r=\x-y\>0, (9) 

X 

where r is the Euclidean distance. Since over¬ 
lap fermions are not ultralocal, the fermion ac¬ 
tion is not strictly reflection positive. Therefore 
one may expect a positive core and a negative tail 

ofC',(r) I22ESI- 

Fig.^shows that there is always a positive core 
near the origin at any value of the cutoff, while the 
tail is only negative for sufficiently large cutoffs. 
Since the topological susceptibility is not affected 
by the cutoff in the number of modes, the growth 
of the positive core and the negative tail must 
compensate each other. 

While the topological charge correlator Cq{x) 
in the continuum should be negative for nonva¬ 
nishing X, the topological susceptibility, which 
is just the integral over the correlator xt = 
f dx (q(O)q(x)} = f dxCq(x) > 0, must be pos¬ 
itive. To solve the clash, formally divergent con- 
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Figure 3. Isosurfaces of topological charge density with \q[x)\/= 0.1, 0.2, 0.3, 0.4, on one timeslice of 
a single 16^ x 32 configuration at /3 = 8.45. The upper (lower) pictures are based on the full (eigenmode 
truncated, with |Im(A)| < 0.076) density. Colour encodes the sign of the surface. 


tact counterterms of the form ci5[x) + C2^8{x) + 
c^A'^Slx) have to be introduced in the continuum 
theory |^ . 

To see whether the expected behaviour is 
reached in the continuum limit, we examine in 
Fig.lHl the a dependence of the charge correlator 
based on the full density. One can clearly see that 
the peak value of the positive core is increasing 
with decreasing a, while the support of the posi¬ 
tive core is shrinking. However, the finite positive 
core at our highest /3 remains at a radius R Ki 2 a, 
indicating the limits of locality of the operator 
q{x) itself. 

In the following we want to further investigate 
the cluster properties of the full density and com¬ 
pare our results obtained on the coarse 12^ x 24 
lattice at /3 = 8.10 and the finer 16^ x 32 lattice at 
(3 = 8.45, which have a similar physical volume. 

In Fig. El (a) we show that the number of 
clusters^ reaches a maximum, which is higher 
nearer to the continuum, at qcut=0.20 (0.25) on 


®The maximal number of clusters of the truncated density 
is an order of magnitude smaller. 


the coarse (fine) lattice. For a cutoff as low as 
9 cut = 0.05 gmax there are, independent of the 
lattice spacing, 2 totally dominating oppositely 
charged clusters I2I|. Fig.El(b) demonstrates how 
the largest cluster approaches a packing fraction 
of 50 %. The behaviour for the second largest 
cluster is similar. 

Fig. El ( c) shows the maximum of the minimal 
distance between each point of the largest cluster 
and every point of the second largest one. The 
small value for low cutoffs in q implies that the 2 
dominating clusters are tangled and intertwined 
in a complicated way. 

In Fig. El(d) we plot the point-to-point cluster 
correlation function fc{R), which measures the 
probability of a pair of points at a distance R to 
be in the same cluster, evaluated at the largest 
possible distance Rmax- We find that the largest 
cluster starts to percolate at g^ut ~ 0.20 gmax for 
the full density independent of the lattice spacing. 

The fractional volume and the total charge of 
all clusters ordered by their size on one typical 
configuration on the 16^ x 32 lattice are plotted in 
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Figure 4. The topological charge correlator on the 12^ x 24 lattice at /? = 8.10. We compare the correlator 
based on the eigenmode truncated version of q{x) for various Acut with the one based on the full density. 
The right figure magnifies the region in R where the eigenmode truncated correlator turns negative. 



Figure 5. The topological charge correlator based on the full charge density on the 12^ x 24 lattice at 
f3 = 8.10 {Nconi = 53) and on the 16^ x 32 lattice at /3 = 8.45 (IVconf = 5). The right figure is a 
magnification of the left to show the negativity of the correlator. 
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Figure 6. Cluster analysis: The gcut dependence of the total number of clusters of full density (a), the 
fractional size of the largest cluster (b), the distance between the 2 leading clusters (c) and the cluster 
correlation function /c(i?max) at the largest possible distance (d) for the full densities. The data is 
plotted both for the 12^ x 24 lattice at /3 = 8.10 and for the 16^ x 32 lattice at /3 = 8.45, which have a 
similar physical volume. Also shown is the fractional volume (e) + (f) and the total charge (g)+(h) of the 
clusters of full charge sorted by their size for gcut = 0.10 q^ax (left) and 0.30 gmax (right) for one single 
configuration on the 16^ x 32 lattice. 










































Figs. ini(e)-(h). Comparing the data for (7cut=0.10 
9 max and 0.30 gmax the striking difference between 
the percolating and nonpercolating regime can be 
observed. 

A preliminary dimensional analysis indicates 
that we are rediscovering for the full density 
the picture of 2 dominating approximately 3D 
sign coherent clusters m at low cutoffs, whereas 
we see ID highly charged small clusters at high 
thresholds. This picture means that the QCD 
vacuum model with 4D coherent (anti)instantons 
with a typical instanton radius of 0.3 — 0.4 fm 
is strongly modified if quantum fluctuations of 
all scales are taken into consideration. The 
eigenmode-truncated density is more compatible 
with the conventional multiple-lump picture. 
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